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ABSTRACT 

Current evolutionary models imply that most Cataclysmic Variables (CVs) have /Lb < 2 hours and are Dwarf 
Nova (DN) systems that are quiescent most of the time. Observations of nearby quiescent DN find that the 
UV spectrum is dominated by the hot white dwarf (WD), indicating that it provides a significant fraction of the 
optical light in addition to the quiescent disk and main sequence companion. Hence, identifying a faint, quiescent 
CV in either the field or a globular cluster (GC) from broadband colors depends on our ability to predict the WD 
contribution in quiescence. We are undertaking a theoretical study of the compressional heating of WDs, extending 
down to very low time averaged accretion rates, (M) ~ 10“' 'Mg yr" 1 , which allows us to self-consistently find 
the T e ff of the WD. We demonstrate here that most of the compressional heating occurs in the freshly accreted 
envelope and that the WD core temperature reaches a fixed value on a timescale less than typical evolutionary 
times. Since nuclear burning is unstable at these (M)’ s, we have incorporated the recurrent heating and cooling 
of the WD core throughout the classical novae limit cycle in order to find the T e s-(M) relations. Comparing to 
observations of field DN confirms the (M)-P OI \, relation of disrupted magnetic braking. We also predict broad-band 
colors of a quiescent CV as a function of (M) and companion mass and show that this leads to the identification 
of what may be many CVs in deep HST images of GCs. 

Subject headings: binaries: close—globular clusters: general—novae, cataclysmic variables—white dwarfs 


1. Introduction 

Dwarf Novae (DN) systems contain a white dwarf (WD) ac¬ 
creting matter at time-averaged rates (M) < 1 0~ 9 M & yr from 
a low-mass (< 0.5 Mq typically) stellar companion (see Osaki 
1996 for an overview). At these (M)’s, the accretion disk is 
subject to a thermal instability which causes it to rapidly trans¬ 
fer matter onto the WD (at M 2> (M)) for a few days to a week 
once every month to year. The orbital periods of these binaries 
are usually less than 2 hours (below the period gap), but there 
are also DN above the period gap, > 3 hours (see Shafter 1992). 
The M onto the WD is often low enough between outbursts that 
the UV emission is dominated by the internal luminosity of the 
WD. Indeed recent spectroscopy has resolved the WD’s contri¬ 
bution to the quiescent light and found effective temperatures 
T e ff ~ 10,000-40,000 K (see Sion 1999). 

The measured internal WD luminosity is larger than expected 
from an isolated WD of similar age (« Gyr), indicating that it 
has been heated by accretion (Sion 1985). Compressional heat¬ 
ing (i.e. internal gravitational energy release) appears to be the 
main driver for this re-heating (Sion 1995). Sion’s (1995) esti¬ 
mate for internal gravitational energy release within the WD (of 
mass M and radius R ) was L ss 0.15 GM{M)/R. However, we 
show in §2 that the energy release actually depends on the ther¬ 
mal state of the WD interior and that the dominant energy re¬ 
lease is in the accreted outer envelope, giving L ss 3 kT c (M)/jim p , 
where // « 0.6 is the mean molecular weight of the accreted ma¬ 
terial, T c is the WD core temperature, m p is the baryon mass, 
and k is Boltzmann’s constant. 

The theoretical challenge that we address in §3 is how to cal¬ 
culate T c as a function of ( M ), and thus find 7' c fr. Because of 
unstable nuclear burning and the resulting classical novae cy¬ 


cle, the H/He envelope mass changes with time, allowing the 
core to cool at low accumulated masses and be heated prior to 
unstable ignition. We use nova ignition to determine the maxi¬ 
mum mass of the overlying freshly accreted shell, and find the 
steady-state (i.e. cooling equals heating throughout the classi¬ 
cal novae cycle) core temperature, T c , as a function of ( M) and 
M. 

We compare our calculations to HST/STIS observations and 
infer ( M) on the timescale of 10 6 years. We find that DN above 
the period gap have (M) « KL 9 M Q yr -1 , while those below 
have (M) ss 10 -u) Mq yr -1 , consistent with that expected from 
traditional CV evolution (e.g. Howell, Nelson, & Rappaport 
2001), even those that involve some “hibernation” (Shara et 
al. 1986; Kolb et al. 2001). The result is more surprising if 
the much weaker magnetic braking laws of Andronov, Pinson- 
neault, & Sills (2001) are correct. We also predict the mini¬ 
mum light (My) of CVs in quiescence for a range of (M), WD 
mass, and companion mass. This assists the search for the pre¬ 
dicted large population of CVs with very low mass companions 
(< O.IMq) that are near, or past, the period minimum (Howell, 
Rappaport, & Politano 1997). Observations already show that 
the WD fixes the quiescent colors of these CVs and our calcula¬ 
tions are useful for CV surveys in the field (e.g. 2DF, SDSS, see 
Marsh et al. 2001 and Szkody et al. 2002) and globular clusters. 

2. Compressional Heating: Core-Envelope 
Contrast 

Compressional heating is the energy released by fluid ele¬ 
ments as they are compressed by further accretion. The impor¬ 
tant feature of this heating mechanism is that the heat is released 
in the WD interior , and thus is radiated on a timescale which 
is longer than the time between DN outbursts. Contrast this 
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to the gravitational potential energy released by the infalling 
matter, (GMm p /R per baryon) which is deposited at, or near, 
the photosphere and is rapidly radiated away. Such infall en¬ 
ergy does not get taken into the star because in the upper at¬ 
mosphere (where T -C T c ) the time it takes the fluid to move 
inward is much longer (by at least T c /T) than the time it takes 
for heat to escape. This means infall energy is not important for 
setting the internal thermal state of the WD; it simply has no 
influence there. Also, once accretion has diminished for longer 
than the time to radiate the infall energy away (such as in DN 
quiescence), it is no longer relevant to the observed luminosity. 

An additional energy source in the WD interior is slow nu¬ 
clear burning near the base of the accreted H/He layer. This is 
significant when the accreted layer becomes thick, eventually 
becoming thermally unstable and leading to a Classical Nova. 
The nova energy is assumed to be radiated away in the explo¬ 
sion, but we have found that slow burning before that point con¬ 
tributes an amount of energy to the WD interior comparable to 
the energy released by compression. Our calculations take ac¬ 
count of both compressional heating and slow nuclear burning 
to determine the WD’s thermal state. 

We now sketch how compressional heating is included in our 
stellar model, demonstrate that heating in the envelope dom¬ 
inates that in the core and estimate its magnitude. This is a 
discussion specific to DN with low (Af), and the reader should 
consult Nomoto (1982) for a complete account. The simplest 
estimate of the compressional energy release is the gravitational 
energy liberated as a fluid element moves down in the WD 
gravitational field, g = GM/R 2 . In the non-degenerate outer at¬ 
mosphere, a fluid element moves a distance of order the scale 
height, h = kT/nm p g, in the time it takes to replace it by ac¬ 
cretion, giving L ~ ( M)gh ~ ( M)kT / pm p . This exhibits the 
correct scaling, notably the dependence on /; which is a con¬ 
trasting parameter between the accreted H/He envelope and the 
C/O core. 

To calculate the actual heat release we consider the local heat 
equation 


ds ds _ _ dL 

T— = T—+TvS7s = + e N , 

dt dt dM, 


( 1 ) 


where e,v is the nuclear burning rate, s is the entropy, and v = 
~(M)r/4nr 2 p is the slow downward advection speed from ac¬ 
cretion. The entropy profile is fixed by the temperature gradient 
needed to carry the luminosity outward and thus we simultane¬ 
ously solve equation (1) with the heat transport equation, using 
opacities and conductivities from Iglesias & Rogers (1996) and 
Itoh et al. (1983) to find the thermal structure of the accreted 
envelope and the outer edge of the C/O core. For an analytic 
understanding, we neglect nuclear burning and d/dt and use 
hydrostatic balance to recast equation (1) into 



Entropy decreases inward (i.e. the envelope and core are not 
convective), so this is an outward L. In the non-degenerate en¬ 
velope, where s = kln(T 3 / 2 / p) / pm p , one more approximation 
is necessary to obtain an analytic form. For an atmosphere in 
which L is constant with depth, the envelope satisfies r 8 - 5 oc 
P 2 . Though L is not constant here, we use this to get an esti¬ 
mate, integrating down to the isothermal core and finding L « 
3 kT c (M)/pm p . 

Now consider the degenerate C/O core. For the (Af)’s and 
typical M = 0.6Mq WD of interest here, the entropy is in the 



Fig. 1.— The Hydrogen/Helium envelope and outer core in temperature and 
pressure for three different values of accumulated mass, M- dCC / Mq = 0.5, 
1.5, and 2.5 for M = 0.6M® and (Af) = 1(T 10 Mq yr -1 . The external surface 
luminosity of the WD and the corresponding T e ff is also listed. The part of the 
star off to the right of the figure (P > 10 20 erg cm -3 ) is the isothermal inner 
core. 


liquid ions at T c « 10 7 K. The time it takes to transport heat 
through the interior is ~ 10 7 yr -C M/(M), so the core is isother¬ 
mal and any compression is far from adiabatic. 1 Due to uncer¬ 
tainty from the classical novae cycle, we don’t know whether 
the C/O core is secularly increasing in mass, but if it were, al¬ 
most all of the work of compression goes into increasing the 
electron Fermi energy. The integrated heat release in the core 
would then be L « \5kT c {M) / p,m () (Nomoto 1982) for a 0.6A/q 
C/O WD, where p, « 14 is the ion mean molecular weight. 

Due to the mean molecular weight contrast between the ac¬ 
creted envelope and the core, the energy release in the core is 
about a factor of five smaller than that in the envelope. Thus, for 
a given amount of compression of the star, the entropy drop for 
material in the accreted layer is much larger than for material 
in the core. Despite its comparatively small mass, the accreted 
layer is the main source of compressional heating. 

3. Finding the Equilibrium Core Temperature 

For this initial study, we dropped the time-dependent term in 
equation (1), and presumed that the C/O core mass was con¬ 
stant throughout the classical nova cycle, thus only accounting 
for the compressional heating and e,v in the accreted layer. This 
method improves on that of Iben et al. (1992) by allowing the 
accreted envelope mass to change through the 10 5 year classical 
nova cycle. Early in the cycle, the mass of the accreted layer is 
small, compressional heating is small, and the WD cools. Later 
in the cycle, the accreted layer becomes thick enough that com¬ 
pressional heating along with slow hydrogen burning releases a 
sufficient amount of energy to heat the core. As the WD has a 
large heat capacity, reaching the equilibrium T c where the heat 
exchanged between the envelope and core averages to zero over 
a single classical nova cycle takes « 10 8 years. Since this time 
is shorter than the time over which (Af) changes due to chang¬ 
ing orbital period, we construct such equilibrium accretors for 
a given M and (Af). 

To do this construction, we first fix T c at the outer edge of the 

1 This is in contrast to the rapid accretion rates (M) 10 yr -1 consid- 

ered for more massive Type la progenitors, where the interior undergoes nearly 
adiabatic compression (see Bravo et al. 1996). 
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Fig. 2.— The solid line is the luminosity at the outer edge of the core and 
the dashed is the surface effective temperature of the WD, both as a function of 
accumulated mass up to the classical nova ignition for the equilibrium model 
of Figure 1, M = 0.6Mg, (M) = l(T 10 Mg yr 1 . Positive luminosity is outgoing 
and epochs of core cooling and heating are indicated. 


C/O core at a pressure high enough so that the changing accu¬ 
mulated mass has little direct effect. With a radiative-zero outer 
boundary condition, we integrate our structure equations with 
equation (1) to find the thermal state for an (M) and accreted 
layer mass. See Figure 1 for examples of the resulting T-P 
relations. We then evaluate the luminosity across the chosen 
location (the right edge of the plot in Figure 1) for different ac¬ 
creted layer masses, up to the unstable ignition which is found 
by comparing the T and p at the base of the accreted layer with 
analytic ignition curves (Fujimoto 1982). 

We vary T c to find an equilibrium model, where the “core 
luminosity” (L core ) averages to zero over the classical nova cy¬ 
cle as shown in Figure 2. The quiescent T e ff for the same cycle 
is also shown in Figure 2. At the nova outburst we assume 
that the accreted shell is expelled, and that, due to the rapid¬ 
ity of this event, it does not appreciably heat the WD. The re¬ 
sulting equilibrium core temperatures for ( M) = 1O _1 °M0 yr~' 
are T C /10 6 K = 9,7.5 and 8.5 for M = 0.4,0.6 and l.OM 0 . The 
OAMq star is hotter than the 0.6 Mq star because it has a larger 
maximum accumulated mass that leads to a longer period of 
core heating. For a 0.6M o WD, the core temperatures are T c j 1 0 f ’K = 
4,5.3,12.2 and 18.0for (M) /Af 0 yr 1 = 10- n ,3.2 x 10- n ,4.2 x 
1(T 10 and HT 9 . 

The T eS during the classical nova cycle varies over a rela¬ 
tively narrow range that allows us to compare to observations. 

For field CVs, the large set of STIS observations by Szkody et 
al. (2001) and previous observations (Urban et al. 2000) pro¬ 
vide spectra of quiescent WDs in DN. These measurements 
are made during deep quiescence when the accretion luminos¬ 
ity is negligible and are long enough after the outbursts that 
other emission mechanisms (e.g. Pringle’s (1988) suggestion 
of radiative illumination of the WD) have faded. The observed 
Teff’s thus measure the heat directly from the WD interior. This 
comparison to observations indicates that below the period gap, 

( M) rs 1O” 1O M0 yr -1 and the WD masses are in the range 0.6- 
1.0 Mq. This agrees with the expectation from Kolb & Baraffe 
(1999), who find ( M) rs 5 x \Q~ xi Mq yr -1 at an orbital period 
of 2 hours presuming angular momentum losses from gravita¬ 
tional waves alone. Above the period gap, the 7)n is higher, and 
we estimate ( M) rs 1 0 9 M 0 yr 1 . For a graphical comparison, 
see Townsley & Bildsten (2001). This general agreement with 


data from field WDs in which the internal luminosity is directly 
visible gives us confidence that our calculations can be applied 
to other quiescent DN systems. 

We predict that a 0.6 Mq WD above the gap has T c = 1.8 x 10 7 
K and, if in equilibrium below the gap, T c = 7.5 x 10 6 K. How¬ 
ever, if the WD does not have time to cool as it traverses the gap, 
it will be hotter than our calculation implies. We estimate this 
cooling time from the current WD cooling law (e.g. Chabrier 
et al. 2000), L coo i ~ 10~ 2 Lq(T c / 1.8 x 10 7 K) 2,5 , along with the 
heat capacity of the core, M3kg/Pitn p , giving At rs 0.5 Gyr. 
Since this is comparable to the estimated time spent in the gap 
(Howell et al. 2001), our equilibrium assumption below the gap 
is likely safe. However, note that about 0.2 Gyrs after accretion 
halts, the WD will enter the ZZ Ceti instability strip! 

4. Application to Globular Cluster Populations 

Due to the high frequency of stellar interactions in Globu¬ 
lar Clusters (GCs), an abundant population of CVs is expected 
to be found there, especially at low ( M). CVs in GCs are com¬ 
monly searched for via the presence of hydrogen emission lines 
or X-ray emission (as recent Chandra observations have found; 
Grindlay et al. 2001a, Grindlay et al. 2001b), and this method 
is fruitful. We show that these systems (as well as CVs cross¬ 
ing the period gap or those “hibernating” post-novae, Shara et 
al. 1986) can also be identified by their position in a color- 
magnitude diagram (CMD). By using our theory of the thermal 
state of the WD, it is possible to predict the broadband colors 
of quiescent CVs. 

An excellent example is NGC 6397 (King et al. 1998; Taylor 
et al. 2001). Figure 3 shows a CMD of NGC 6397 with our ini¬ 
tial results. The data points are objects which meet the proper- 



FlG. 3.— Color-Magnitude Diagram of NGC 6397. The dots are plausible 
cluster members that are below the MS (HST observations by King et al. 1998) 
and the 0’s are the “Non-Flickerers” from Taylor et al. (2001). The MS line is 
from Baraffe et al. (1997) for the cluster [M/H] of-1.5 and an age of 10 Gyr. 
The dashed line is from Bergeron et al. (1995) for DA WDs with logg = 8. The 
lines connecting the MS to the WD sequence are our current calculations of 
the WD+MS binary at the specified ( M). The highest 70 during the classical 
nova cycle has been used for the WD in each case (see Figure 2). No disk has 
been included. All curves have been put at the distance and reddening of the 
cluster, (m-M), = 12.05 and E(V-I) = 0.288. 
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motion criteria for cluster membership and which are below the 
MS. The lines were produced by superposing a WD with the 
maximum 7 df for the indicated (M) with a MS star. Due to un¬ 
certainty in the theory of quiescent disks (Menou 2000), no disk 
contribution has been added. Note, however, that a constant 
T ~ 5000 K disk like those indicated in eclipse maps (Wood 
& Crawford 1986) would have a V — I color of 1.24 includ¬ 
ing the cluster reddening. Except for near the WD cooling line 
(dashed curve), where the WD dominates, the 7 magnitude is set 
by the MS companion. The large dots along the 10 9 M 0 yr 1 
and 10“ 10 Mq yr lines indicate where the MS companion is 
0.3, 0.2, 0.15 and 0.1 Mg, and two additional points at 0.09 and 
0.085 Mq are indicated on the 10 ~ n M & yr~' line. This imme¬ 
diately provides a number of candidate systems (namely, data 
in this part of the CMD). 

The circled points are the “non-flickerers” (Cool et al. 1998) 
recently reported by Taylor et al. (2001). The three at 7« 22.25 
are very strong Ha absorbers (consistent with a DA WD) and 
were not detected by Chandra (Grindlay et al. 2001b). These 
authors had discussed these systems as possible helium WDs 
with millisecond pulsar companions, though, given our work, 
we would claim that these are hot WDs with « 0.15M 0 MS 
companions. In addition, the population of data points in this 
diagram with respect to our theoretical curves will eventually 
constrain CV evolutionary scenarios. If we assume many of the 
data points are CVs, we already see that most systems with high 
(M) have O.15-O.3M 0 companions. If confirmed as members 
of the cluster, the stars below the (M) = 1CT 1 'Af@ yr -1 line could 
well be the long-sought post-turnaround systems with (M) = 
10 -12 Mq yr -1 and companion masses < 0.09M 0 (Howell et al. 
1997). 

5. Conclusions and Discussion 


envelope is thin (i.e. right after the classical nova) is balanced 
by that regained when the envelope is thick. This equilibrium 
determines the 7(n of the WD throughout the classical nova cy¬ 
cle. Our models agree with the observations of dwarf novae in 
deep quiescence and imply (M) « 1CT 1O M 0 yr -1 just below the 
period gap and (M) ss 1CT 9 M 0 yr~* just above the period gap 
for WD masses in the range 0.6-1.0 Mq. 

Our 7)if calculations provide a prediction of the colors of qui¬ 
escent DN. Using MS stellar models, we have predicted where 
a DN should appear in a color-magnitude diagram as a function 
of (M) and the mass of its companion. Many unidentified ob¬ 
jects appear in the relevant regions of the detailed CMDs which 
have been obtained for globular clusters by HST. The number 
of such systems in the field will increase due to upcoming sur¬ 
veys (such as SDSS and 2DF, see Marsh et al. 2001), and will 
push to lower ( M) systems. 

Though our initial efforts have met with apparent success, 
there is still much to be done. We need to vary the metallic - 
ity of the accreted material, lowering to values appropriate for 
globular cluster science. This could change our results at large 
( M), but at low (M)’s, the ignition mass is set by pp burning 
and will likely not change too much. We also need to relax our 
initial assumptions, e.g. by including WD excavation or accre¬ 
tion and accounting for thermal evolution of the WD. 

The internal thermal state of the WD has been a longstand¬ 
ing uncertainty in classical nova work, as has the question of 
how much mass is ejected in the explosion (Gehrz et al. 1998). 
Our work provides the first calculation of the internal thermal 
state of a WD undergoing classical novae, and will eventually 
lead to self consistent calculations for ignition masses, includ¬ 
ing variations of the metallicity. This will be an improvement 
on previous work (e.g. Prialnik & Kovetz 1995) which treated 
1\ and ( M) as two independent parameters. 


We have evaluated the action of compressional heating on 
accreting WD interiors and shown that most of the compres¬ 
sional energy release takes place in the accreted envelope, and 
is thermally communicated to the core. The maximum envelope 
mass is set by the unstable nuclear burning that causes a clas¬ 
sical nova runaway and most likely expels the accreted mass. 
We have constructed equilibrium accretors which have constant 
core temperatures such that the heat lost from the core when the 
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